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Abstract 

We show that the contact process in an ensemble with conserved total 
particle number, as simulated recently by Tome and de Oliveira [Phys. 
Rev. Lett. 86 (2001) 5463], is equivalent to the ordinary contact pro- 
cess, in agreement with what the authors assumed and believed. Sim- 
ilar conserved ensembles and equivalence proofs are easily constructed 
for other models. 
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A statistical mechanical ensemble is a probability distribution on a set 
of microstates. For systems in thermodynamic equilibrium the ensemble 
is known a priori. The equivalence of the standard equilibrium ensembles 
(microcanonical, canonical, ... ) is a well-known fact. In Monte Carlo sim- 
ulations an easy way to generate a given ensemble is to impose detailed 
balancing on the transition rates between the microstates. A system of 
statistical physics may be defined more generally, and more implicitly, not 
by an a priori given ensemble but by the transition rates between the mi- 
crostates. The stationary (= time independent) solution of the resulting 
stochastic equation of motion then constitutes an ensemble which, unlike 
the standard equilibrium ensembles, cannot usually be expressed in explicit 
analytic form. Such stationary ensembles Q arise, e.g., in reaction-diffusion 
type systems, of which one example is the object of study here. 

We consider the contact process (CP) introduced by Harris which 
takes place on an infinite regular lattice each of whose sites may be empty 
or occupied by a particle. The transition rates are such that particles (i) 
are created on an empty site at a rate n^ri/z, where nnn is the number 
of occupied nearest neighbors of that site and z the lattice coordination 
number; and (ii) annihilate spontaneously and independently at a rate k. 

The contact process is of interest because its stationary ensemble exhibits 
a phase transition between an active state (particles present) and an absorb- 
ing state (empty lattice) for a certain critical value k = kc of the annihilation 
rate [^. Monte Carlo simulations of the active state on finite lattices and 
for k close to kc suffer from the complication that in any long enough run 
there will inevitably occur a fluctuation that pushes the system irreversibly 
into the absorbing state. Various procedures have been proposed for dealing 
with this "accidental death" , which occurs similarly in other model systems 
with transitions to absorbing states. Very recently Tome and de Oliveira 
[Ql had recourse to an elegant method, used earlier in a different context by 
Ziff and Brosilow Q, and which circumvents the accidental death problem 
altogether. Tome and de Oliveira replace the ordinary contact process (CP) 
by a version with strictly conserved particle number N. In this conserved 
contact process (CCP) particles (i) are created at the same rate as in the 
CP; but (ii) every creation is accompanied by the simultaneous annihilation 
of another particle picked at random on the lattice. It is our purpose here to 
show that the conserved ensemble (CCP) is equivalent to the nonconserved 
ensemble (CP). This equivalence was assumed and believed to be true by 
the authors Q, who nevertheless feel the need to invoke their results as an 
additional argument. The proof is simple and, just like the conserved en- 
semble method itself, can easily be adapted to other problems. 

Let j denote a lattice site and let rij be its occupation number {nj = or 
= 1). It will be convenient to work with "spins" sj = 2nj — 1, so that sj = 
±1. The symbol s = {sj} then represents the set of occupation numbers. 
We associate |^] with each s a basis vector \s) and with the time dependent 
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probability distribution P{s,t) the state vector 

\p)t = Y.ns,t)\s) (1) 

s 

The master equation for P{s,t) is then equivalent to the evolution equation 

d\P)t 



dt 



for the state vector, in which the "master operator" W is the infinitesimal 
generator of the transitions. The Sj = ±1 may be considered as the eigen- 
values of a Pauli spin matrix (t|. Since <tJ reverses sj, one may express W 
in terms of the crj and (t| |^]. For the standard CP with annihilation rate k 
this yields W = with 

= Were + kW^nn (3) 

Wee = EE'K-D^^i^ (4) 

i & 

where j + 5 is a neighbor site of j and denotes 1/z times the sum 
on all neighboring sites. For the CCP with exactly N particles one finds 
W = vyccp with 

w'^^^ = ^EE'EK< - (6) 
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(in which the term i = j is identically zero) . The question is now to demon- 
strate the equivalence, in a sense to be appropriately defined, of W^^ and 
yyccp_ Q^j. pj-QQf proceeds by two steps. The first one is to rewrite expres- 
sion (^) for the operator W*^^^ by means of the substitution 

ctJ< - 1 = {a] - 1) + « - 1) + (aj - 1)« - 1) (7) 

Using that in the conserved ensemble = J2ii^ + )/2 we find 

W°^^ = Were + Wann^ + Wr, (8) 

Here k = TVio /N in which 

is the operator for the total number of pairs of neighboring sites of which 
one is occupied and the other empty; and the "remainder" Wr, is given by 

^ J- 5 i III ^^Q^ 
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We remark that whereas Were and Wann are themselves master operators 
(i.e. correspond to master equations), the remainder Wr is not; the reason 
is that although it conserves the total probability, it does not conserve the 
positivity of an initial probability distribution. 

The second step of the proof is to consider the time evolution of the 
ensemble average of an arbitrary product / = af_^ af^ . . . af^, knowing that 
all physical observables are linear combinations of such products. Let {f)^^^ 
denote the average of / in the conserved ensemble. It may be calculated as 
the scalar product 

(/)ccp = ^o\f\P)t (11) 

where \P)t = exp(W°'^^t)|P)o and where the "projection state" (0| is de- 
fined by (0| = X]s('5| with the sum running on all s irrespective of the 
number of particles present. From Eqs. ([TT|) and (P) it follows that 



J / f\CCP 

= (/Werc)^" + (/Wann/^)^^ + (/Wr),°^^ (12) 

In the usual nonconserved ensemble one has in the same way from Eqs. (|ll| ) 
and (|3|) that 

We now compare Eqs. (^) and ([I^). In the limit of an infinite system (that 
is statistically invariant under translations) the operator k will have van- 
ishingly small fluctuations around its average (k)^*^^ = (A/io)f'^^/-/V. Hence 
we may replace k by this average. Next we will show that the last term 



on the RHS of Eq. (12) vanishes in the limit of an infinitely large system. 
It is not sufficient for that to invoke the factor 1/N in expression ( [ToD for 
Wr. Instead, we will make the average (/Wr)^*^^ fully explicit. In order to 
do so we commute the crj through / using that [cr|,(Tj]^ = 0, and use that 
(0|a^ = (01. This gives 



^ p=lq=l 5 ' " " " ■ ' (14) 

The RHS of this expression involves a sum on terms and the summand 
makes no reference to the system size. Therefore it is now clear that in view 
of the prefactor 1/A^ the average (^) vanishes in the limit of infinite system 
size. Hence Eq. ([T^ ) may be replaced with 

J / f\CCP 

= (/Were)^" + (/Wann)r (15) 



Comparison of (13) and (^) shows that any spin average evolves accord- 
ing to the same equations of motion in both ensembles provided we have 
(k)^*^^ = fc, i.e. 

{Ni,)rvN = k (16) 
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This equation relates the CP parameter k to a time dependent CCP average. 
It can be satisfied only if that average is time independent, hence equal to 
its stationary state value. Therefore the condition for (|l^) and (12) to be 
equivalent becomes 

(-A/'io)^t"ar/A^ = k (17) 

where the average is on the stationary CCP ensemble. In the stationary 
state of the CCP, or for fluctuations around it that in the large-A^ limit 
affect (A/io)stat /-^ negligibly, equality (p!7|) guarantees that the two averages 
{f)t^ and obey the same equations. It follows in particular that for 

any / we have (/)^tat — (/)st'at ■ This establishes the equivalence of the two 
stationary ensembles. A notable case which on the basis of the preceding 
discussion is excluded from the equivalence, is the relaxation towards equi- 
librium of an initial state in which all particles occupy random positions: 
In such a nonequilibrium process ffl the average of TVio is not constant (in 



either ensemble) and (17) cannot be satisfied. 



The equivalence condition ([T7|) derived here was introduced by Tome and 
de Oliveira (Eq. (3) of Q), who exploit it to determine the critical behavior 
of the contact process. Here we have furnished the proof (albeit a physicist's 
one) that their procedure was correct. 

The construction of conserved ensembles by simultaneous execution of 
elementary transitions at uncorrelated lattice locations can easily be adapted 
to other models of interest in statistical physics. The same is true for the 
equivalence proof of this note. 
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